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ABSTRACT 

The  objective  of  this  investigation  is  to  derive  the  reliability  and  the 
associated  probabilistic  failure  criterion  for  composite  materials  under 
combined  stress.  In  the  analytical  derivation,  the  concept  of  joint  probability 
was  used  and  applied  to  the  Weibull  distribution  function.  In  applications,  given 
the  experimental  measurements  of  the  necessary  statistical  parameters  for  the 
specific  composite,  the  probabilistic  criterion  of  the  composite  failure  and  the 
reliability  of  the  specific  structure  can  be  predicted. 

Graphical  representations  for  the  joint  reliability  and  joint  failure  contours 
were  made  in  two  and  three  dimensional  space  for  the  several  different  sets  of 
statistical  strength  parameters  to  illustrate  the  effect  of  parameters  on  reliability. 
Such  understanding  will  enhance  selection  of  fiber  and  matrix  (which  have  their 
own  statistical  strength  parameters)  and  can  lead  to  improvement  in  reliability  of 
some  composite  components  in  an  aircraft  structure.  These  reliability  and 
failure  concepts  can  also  be  used  in  repair  problems  by  selecting  the  proper 
composite  material  with  the  appropriate  statistical  parameters. 
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I.   INTRODUCTION 

A.   FAILURE  AND  RELIABILITY  OF  THE  COMPOSITE  MATERIAL 

In  a  composite,  the  matrix  and  reinforcement  materials  are  combined  to  form 
a  macroscopically  multiphase  material  .  The  macroscopic  state  gives  rise  to  a 
material  which  is  internally  redundant  against  local  defects  and  damages. 
[Ref.1]  The  ply  is  the  basic  unit  in  constructing  composites  which  consists  of 
fibers  embedded  in  a  matrix.  It  is  possible,  therefore,  to  relate  composite  design 
properties  to  corresponding  constituent  properties.  As  a  result,  it  is  possible  to 
design  structural  components  made  of  fiber-composite  when  the  design 
requirements  and  the  properties  of  the  candidate  constituent  materials  are 
known.  For  high  performance  structures,  the  unidirectional  fiber-matrix 
composites  consist  of  strong  stiff  fibers  embedded  in  a  comparatively  ductile 
matrix  .  For  a  rational  structural  design,  the  strength  of  such  unidirectional  ply 
needs  to  be  expressed  in  terms  of  the  stress  tensor.  Because  the  ply  is 
composed  of  numerous  fibers,  the  composite  failure  process  is  sequential, 
therefore  ad  hoc  failure  criteria  can  not  be  applied  to  analyze  the  composite 
failure.  Furthermore,  such  conventional  failure  criteria  such  as  Mises  and 
Tresca's  criteria  are  related  only  to  a  single  parameter  (mean  function  of 
failure).  As  a  result,  they  can  not  be  applied  to  multiple  parameter  functions,  so 
suitable  general  relations  to  represent  failure  contours  of  composite  materials 
are  required. 

In  the  application  of  a  one-dimensional  failure  criterion,  the  magnitude  of 
stress  (one  component)  at  each  spatial  location  of  the  structure  is  mapped  into 
respective  points  on  the  stress  space  (which  is  a  line).  If  all  the  points  within  the 
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domain  are  bounded  by  the  tensile  strength  limit  and  the  compressive  strength 
limit,  then  the  entire  structure  may  be  considered  safe.  The  reliability  under 
unidirectional  stress,  then,  is  the  product  of  the  reliability  of  each  spatial 
segment  at  the  respective  stress  level.  Under  planar  combined  stress,  the 
magnitude  of  stresses  (3  components)  at  each  spatial  location  of  the  structure  is 
mapped  into  respective  points  on  the  six  dimensional  stress  hyperspace.  So  the 
stress  that  is  interior  to  a  domain  bounded  by  the  failure  surface  is  safe  and  the 
failure  criterion  is  the  analytical,  graphical,  or  numerical  representation  of  the 
failure  surface.  To  formulate  the  failure  criterion,  consistency  of  mathematical 
operations  can  be  assured  by  adhering  to  the  quotient  rule  of  tensor  variables. 
Historically,  the  formulation  of  many  failure  criteria  were  intuitively  taken  to  be 
based  on  failure  mechanisms  but  in  fact,  they  are  phenomenological  because 
of  the  inconsistency  between  macro  and  microscopic  failure  modes. 
Phenomenological  failure  criterion  may  be  considered  as  a  mathematical 
model  of  the  material  transfer  function  relating  the  external  excitation  (stresses) 
to  the  material's  response  (failure).  The  phenomenological  mathematical 
model  is  intended  to  aid  experimental  design  to  facilitate  interpolation, 
correlation,  and  retrieval  of  experimental  observations.  They  do  not,  in  general, 
address  the  statistical  variability  of  strength. 

For  a  combined  stress  failure  mechanism,  when  one  stress  component  does 
not  affect  the  strength  of  another  component,  the  failure  is  mechanistically 
uncoupled;  conversely  when  one  stress  component  affects  the  strength  of 
another  component,  the  failure  is  mechanistically  coupled.  In  addition,  if 
probabilistic  failure  is  also  taken  into  consideration,  the  condition  for  the  failure 
mechanism  is  not  fixed;  it  depends  on  the  stochastic  combinations  of  the 
intrinsic  strengths.     Under  combined  stress,  mechanistic  uncoupling  with 


statistic  independence  can  interact  to  produce  phenomenologically  coupled 
contours.  On  the  other  hand,  mechanistically  coupling  with  statistical 
dependency  can  also  produce  phenomenologically  coupled  statistical  contours. 
These  cases  can  only  be  distinguished  by  examination  of  both  the  contour  in 
the  stress  space  and  the  statistical  strength  distribution  in  the  probability  space. 
In  any  phenomenological  failure  criterion,  the  shape  of  the  failure  envelope 
is  never  completely  known  until  experiments  are  performed  for  all  possible 
combined  states  of  stress.  As  a  result,  complete  experimental  determination  in 
two  or  more  dimensions  place  an  impractical  demand  on  time  and  resources 
because  of  the  large  number  of  tests  needed.  However,  with  the  establishment 
of  a  probabilistic  failure  criterion,  this  experimental  failure  contour  can  be 
related  to  the  shape  parameters  (the  variability)  in  the  material's  principal 
direction.  With  such  supplemental  information,  failure  contours  connecting  the 
mean  value  can  be  established  by  a  smaller  number  of  samples.  Therefore  an 
analytical  model  can  complement  the  failure  characterization  by  reducing  the 
required  number  of  experimental  measurements,  providing  bases  for 
comparison  of  the  material's  performance,  and  simplifying  data  and  data  base 
reduction.  In  structural  analysis  applications,  the  reliability  of  a  specific  structure 
under  a  combined  stress  can  be  determined  by  mapping  the  spatial  domain  into 
the  stress  domain  through  stress  analysis.  Using  the  probability  failure  criterion, 
the  reliability  of  each  spatial  location  corresponding  to  the  respective  stress 
tensor  and  magnitude  can  be  calculated.  The  joint  reliability  of  the  individual 
spatial  location  can  in  turn  be  combined  to  predict  the  reliability  of  the  entire 
structure. 


B.  OBJECTIVE  OVERVIEW 

The  objective  of  this  investigation  is  to  derive  an  analytical  model  for 
composites  failure  and  reliability  under  combined  stress  conditions.  The 
Weibull  distribution  function  was  used  in  the  probabilistic  modeling.  Chapter  2 
presents  the  theoretical  formulation  for  failure  and  reliability,  and  the  resulting 
equations  are  shown  graphically  in  the  following  chapters.  Chapter  3  shows 
the  joint  reliability  and  linearized  failure  Cumulative  Distribution  Function  (CDF) 
in  two  and  three  dimensional  space.  The  effects  of  parameters  which  depend 
on  the  material  were  analyzed.  Chapter  4  shows  the  joint  failure  CDF  in  two 
and  three  dimensional  space.  The  effects  of  material  parameters  were  also 
analyzed.  In  Appendix  B,  these  probability  distributions  were  compared  to 
experimental  data  availabe  in  the  literature. 

Finally,  Chapter  5  provides  conclusions  and  some  remarks  for  the  future 
work  in  this  area. 


II.  THEORETICAL  FORMULATION 

A.  RANDOM  VARIABLES 

A  failure  criterion  for  a  solid  prescribes  the  risk  of  failure  in  terms  of  the 
causes  of  the  failure.  Mechanical  catastrophic  rupture  is  considered  herein. 
The  effect  of  mechanical  failure  may  be  considered  to  be  caused  by  either 
stress  or  stain  depending  on  the  microscopic  failure  processes.  If  the 
constitutive  relation  of  the  material  is  available,  conversion  between  stress  and 
strain  can  be  readily  made  and  the  two  variables  may  be  considered  as 
equivalent.  Similarly,  a  probabilistic  failure  criterion  which  describes  the 
probable  risk  of  failure  can  be  expressed  in  terms  of  stress  or  strain  as  random 
variables.  Stress  is  selected  as  the  random  variable  herein  with  the 
understanding  that  the  variable  can  be  changed  to  strain  with  the  appropriate 
constitutive  relation.  We  denote  the  probabilistic  random  variable  as  Xj,  for  the 
scalar  component  of  the  stress  tensor  ,  and  xj  for  the  realized  random  variable 
(i.e.  the  strength  Xj  under  the  respective  stress  component  oj). 

B.  PROBABILISTIC  CONDITIONS  OF  FAILURE 

Under  combined  stress  conditions,  when  more  than  one  component  of  the 
stress  tensor  assumes  non-zero  values,  the  effect  of  each  component  of  the 
stress  tensor  contributes  to  the  failure  of  the  composite.  Failure  of  the 
composite  does  not  occur  if  failure  by  each  individual  stress  component  does 
not  occur. 

For  the  composite  subjected  to  a  planar  state  of  stress,  we  denote  the  failure 
events  as: 


A:   Failure  caused  by  normal  stress  along  the  fiber  direction  -  a-| 
B:  Failure  caused  by  normal  stress  perpendicular  to  the  fiber  direction  -  o>> 
C:  Failure  caused  by  shear  stress  along  the  fiber  direction  -  og. 
Then  the  failure  of  the  composites  is  either  by  A  or  by  B  or  by  C  or  by  all  three  as 
illustrated  in  Figure  2-1: 

P(A  +  B  +  C)  =  P(A)  +  P(B)  +  P(C)  -  {P(AB)  +  P(BC)  +  P(CA)}  +  P(ABC) 

(2-1) 
As  a  special  case,  if  all  three  events  are  mutually  exclusive,  then  Equation 
(2-1)  can  be  simplified  by 

P(A  +  B  +  C)  =  P(A)  +  P(B)  +  P(C)  (2-2) 


Figure  2-1  :  Joint  Probability 


C.  JOINT  DISTRIBUTION  FOR  FAILURE 
1.   Discrete  Case 

In  the  development  of  the  analytical  expressions  for  the  probability  of 
failure,  we  first  developed  the  general  form  of  the  distribution  for  three  random 
variables  using  the  notation  Xj  for  the  random  variables  and  xj  for  the  realized 
random  variables.  For  the  planar  stress,  the  customary  notation  for  the  range  of 
i  =  1,2,6  is  used.  These  general  expressions  are  then  specialized  for  the 
combined  stress  cases  by  replacing  the  realized  random  variables  by  Oj. 

The  probability  of  failure  under  the  combined  influence  of  three  random 
variables  can  be  expressed  by  the  joint  probability  density  function  (pdf)  of  the 
individual  realized  random  variables: 

P(X,  =  a1,X2  =  a2,X6  =  c6)=f(o1,o2,a6)  (2  3) 

where  f(a-| ,  02,  ag)  >  0  and  ZSZ  f(a-j ,  02,  og)  =  1 .  If  random  variable  Xg  is  to 
take  on  any  one  of  the  values  of  051 ,  og2>  °"63- ••'  °6n  -  ^en  ^e  Jomt  P^f  at 
some  combination  of  X-j  =a-|  j,  X2=°"2j'  and  Xg=agk  can  be  expanded  by 

P(X1  =  a1.,X2=a2.,X6  =  a6k)  =  f(a1.,a2j)a6k)  ^ 

n 
P(X1  =  a1.,X2  =  a2A=f(a1.,a2)  =  If(a1j)a2,a6k) 

'  /      k=o  (2-5) 


1 

P(X2=02j'X6  =  °6k)=f(a2ra6k)  =  ^f(aii'a2ra6k) 


i=0 

m 


(2-6) 


P(X1  =  01i'X6  =  a6k)=f(aii'°6k)=^f(aii'a2ra6k) 

J=o  (2-7) 

and  the  marginal  probability  function  can  be  obtained  for  the  three  random 
variable  case: 


m     n 


P(X1  =  a1i)  =  f(a1.)  =  SIf(a1i,a     aek) 

j=ok=o  (2-8) 


I      n 


P(X2  =  o2j)  =  f(a2j)  =  IXf(cli,c2j,a6k) 


i=Ok=0 

I     m 


(2-9) 


P(X6  =  °6k)  =f  (a6k)  =  SXf(a1j,a2j,a6k) 

i=Oj=0  (2.10) 

which  can  be  written 

If1(a1i)  =  1,if2(a2.)=1)if6(a6k)  =  1 

i=o  j=o  k=o  (2-11) 

I     m     n 

IXIf^.c     a6|\=1 

i=oj=ok=o  (2-12) 

2.   Continuous  Case 

If  three  random  variables  are  considered,  then  Equation  (2-4)  through 
(2-7)  can  be  expanded  by 

P(X,<a1,X2<a2,X6<a6)  =  F(a1,a2,06) 


a      a      o 

6  r      2  r      1 


=  [  [    f  (u,v,w)dudvdw 

(2-13) 


•  oo      —  oo      —  oo 


a      a 

2  r     1 


P(X)<cl,X2<a2)=F(o1,a2)  =  J     J    Jf(u,v,w)dwdudv 


-OO       — oo       — oo 


o      a 

6  r      2 


P(X2<02,X6<a6)=F(o2,a6)  =  J     J     Jf  (u,v,w)dudvdw 

—  oo       — oo       —  oo  \^—~   '  ^/ 

a      0      ~, 

c  1 

P(X1  <  arX6  <  o6)  =  F(<Jra6)  =  J     J    Jf  (u,v,w)dvdudw 


—  oo      — oo      — oo 


(2-16) 
and  the  marginal  distribution  can  be  obtained  by 

°  oo  oo 

P(X,<  a,)  =  F(CTl)  =  \  '/    JfCu.v.wJdvdwdu 


—  OO       — OO       —  oo 
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P(X2  2  a2)  =  F(o2)  =  J  'J    Jf(u,v,w)dwdudv 


a 


P(X6£06)=F(°6)  =  /       J  _  Jf_(U.V,W)dudvdw 


—  oo       — oo      — oo 


which  can  be  written 

oo         oo  oo 

Jf(u)v,w)dudvdw=  1 


(2-20) 

From  the  definition  of  conditional  probability,  the  conditional  probability  of 
02  given  that  o^  occurs  can  be  represented  by  a  function  form: 

U    V  (2-21) 

where  f(o-| ,  02)  =  P(X-|  =o-| ,  X2=c>2) 
a.   General  Case 

If  we  consider  the  combined  stress  case,  the  random  variables  X-|, 
X2,  and  Xq  can  be  assumed  to  be  the  failure  strength  in  the  1 ,  2,  and  6  direction 
respectively.  And  under  a  combined  stress  condition,  the  failure  occurs  when 
either  failure  strength  X-| ,  X2,  or  Xq  is  reached.  Therefore  the  probability  of 
failure  under  combined  stress  X-|  and  X2  can  be  represented  by 

P(X1+X2)  =  P(X1)  +  P(X2)-P(X1,X2) 
and  then  the  above  relation  can  be  related  to  the  joint  distribution  function  as 
follows 

F(o1 ,  o2)  =  P(X1  <  Oi ,  X2<  o2)  =  P(Xi  +  X2;  Xi  <  ^ ,  X2  <  o2)  (2-22) 
and  similarly  for  three  random  variables  case: 

F(a1,a2,a6)  =  P(Xi  <cj1  ,X2<a2,X6<a6)  =P(X1+X2  +  X6  : 
X-j  <a1  ,X2<a2.x6^a6)  (2"23) 


Equation  (2-22)  is  proven  in  Appendix  A-1  using  the  Weibull 
distribution  function  for  the  independent  case. 

If  the  joint  distribution  of  failure  under  the  combined  stresses  of  g-| ,  o>>, 
and  gq  is  considered,  the  strength  failure  criterion  of  a  composite  is  the 
minimum  of  the  failure  strength  X-j ,  X2,  and  Xg,  that  is,  at  Xc  =  X-|AX2AXg.  (  We 
are  not  dealing  with  the  minimum  magnitude  among  X-|,  X2,  and  Xg  but  rather 
one  value  among  X-j ,  X2,  and  Xq  which  is  less  than  the  intrinsic  failure 
strengthes  of  o-j ,  02,  and  og  becomes  the  minimum  value  here.)  So 

Havc2,ae)  =  P{Xc=X1AX2AX6over(a1,a2,a6U 

where  (  o-| ,  02,  gq  )  is  the  three  dimensional  stress  domain.  If  the  relations 
among  c-|  ,  o>>  and  gq  are  given  by 

a1  =B12o-2.  ^2  =  B26^6-  a6  =  B61°"1  (2-24) 

then,  the  joint  pdf  P(X-]  <  a-|  ,X2<G2,Xg<Gg)  can  be  evaluated.  If  o*g  is  assumed 
to  be  zero,  then  the  joint  probability  functions,  P(X-|<g-|,X2<G2)  can  be  divided 
into  two  probability  distribution  functions  corresponding  to  Si  and  S2  domains 
as  shown  in  Figure  2-2.  The  domain  S-|  is  applied  for  the  case  when 
[X-j  /  X2  <  B-|  2  anci  X-j  <  g-|]  which  the  composite  fails  by  failure  strenth 
X2  =  X«|AX2,  whereas  the  domain  S2  is  applied  for  the  case  when 
[X1/X2  >  B-|  2  anci  X2  ^  ^2]  which  the  composite  fails  by  failure  strength 
X-|  =  X-|AX2-   So  the  following  joint  relation  can  be  obtained: 

P(X-|  <g-|,  X2<g2)  =  P(XVX2<Bi2andXi  <o-j) 

+  P(X-|/X2  >  B12  and  X2  <  c2)  (2-25) 

Applying  Equation  (2-25)  to  domains  described  in  Figure  2-2,  then 

P(X1<g1,X2<o2)  =  J   Jf(u,v)dudv  +  J    Jf(u,v)dudv 


s;  -s2 
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Using  end  limits  and  considering    P(X-|<g-|  ^^o^)  =  F(a-j  ,a2),then  the  joint 
distribution  function  for  the  general  case  can  be  represented  by 


F(oi;a2)  =  J    J^_f(u,v)dvdu  +  J     J      f(u,v)dudv 


0         B 


12 


where  u  and  v  are  dummy  variables. 


B    v 

12 


(2-26) 


(T2 


Figure  2-2  :  Domains  For  The  Joint  Distribution 

b.  Specific  Case 

The  joint  probability  of  the  bi-axial  case  (Equation  (2-25))  can  be 
expanded  for  specific  cases. 

(1)  Independent  Case.  If  the  events  X-|  =  o~i  and  X2  =  02  are 
independent  for  all  o-|  and  02,  then  the  joint  probability  can  be  represented  by 
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P(Xi=o1fX2=o2)  =  P(X1-o1)  P(X2=a2) 

or  equivalently 

f(a1,a2)  =  f1(a1)f2(a2)  (2-27) 

If  three  independent  random  variables  are  considered: 
P(X1=sa1,X2=a2,X6=a6)  =  Ppc^)  P(X2=o2)  P(X6=o6) 
f(o1.<^.a6)-f(o1)f(02)f(a6)  (2-28) 

Similarly,  if  the  random  variable  X-|  and  X2  are  independent  for  all 
a-\  and  o2: 

P(Xi<a1,X2<o2)  =  P(X1<o1)  P(X2<a2) 
or  equivalently: 

F(a1,o2)  =  F1(a1)F2(a2)  (2-29) 

where  F-|  (ct-j  )  and  F2(o2)  are  the  marginal  distribution  functions  of  X-|  and  X2 
respectively.  Similarly,  for  the  three  independent  random  variable  case: 

F(X1<o1>X2<a2,X6<a6)  =  P(X1<o1)  P(X2<a2)  P(X6<a6) 
or  equivalently 

F(CT1  ,g2,c6)  =  F(01 )  F(c2)  F(cr6)  (2-30) 

So  the  joint  distribution  function  can  be  expanded  a  step  further  for 
the  independent  case. 

From  Equation  (2-26)  (refer  to  Appendix  A-2) 
F(a1,a2)  =  Fl(ol)  +  F2(a2) 


a  o 

1  ^     i,      n  r      2 


"/     f1(U)F2^]dU-l       f2<V)Fl(B12V)dV 

0  VD12J  o  (2-31a) 

Equation  (2-31  a)  can  be  further  expanded  by  integration  by  parts 

F(aV°2)  =  F,(°l)  +  F2(°2)  "  2Fl(°l)F2(°2) 

12  o  \D^2)  °  (2-31  b) 
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Comparing   Equations   (2-31  a)   and   (2-31  b)   reveals  that  the 
arguments  of  the  pdf  and  CDF  are  exchanged  so  either  Equation  (2-31  a)  or 
(2-31  b)  can  be  used  depending  on  simplicity  of  algebra. 

(a)  Weibull  Distribution  Model.   If  the  Weibull  function  is  applied  to 
the  independent  case,  then  Equation  (2-31  a)  can  be  expanded  by 


F(°1>°2) 


1  -  exp* 


a 


'V 


[  101 


+ 


1  -  exp* 


r°rVku_ 

~JoUAPi 


«r1 


exp 


°0 


1  -  exp 


V 


a    -1 


^«oVu 


~'o    Ihkh 


expi  - 


r  u 


a   -> 
2 


1  -  exp 


B12P; 


rB12u^ 


P. 


du 


du 


V    H    J    )S 

By  rearranging  terms,  the  joint  distribution  function  for  the 
Weibull  model  can  be  obtained  by 


■^•■Kftfr 


«r1 


«. 


B 


u 


f      U       A 
lB12P2J 


V 


du 


<mt 


a    -1 

2 


a 


exp 


P. 


du 


(2-32) 


(2)  Independent  And  Identical  Case.    Equation  (2-31  a)  can  be  further 
expanded  if  the  identical  case  is  considered. 

For  the  identical  case,  F(a)  =  F-|  (a-| )  =  F2(o>>)  ancl  equivalents  f(o) 
=  f-|  (a-| )  =f2(o"2)-  so  expanding  Equation  (2-31  a) 

c  a 

^1.a2)  =  F(a1)  +  F(a2)-J    f  (u)ff^-\iu  +  J     fMF(B,2v)dv 


12 


(2-33a) 


and  similarly  Equation  (2-31  b)  can  be  represented  by 
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F(a1  ,a2)  =  F(a1 )  +  F(o2)  -  2  F^ )  F(o2) 


+ 


fcX^M(^+B"J."W(B,«¥)dw 


(2-33b) 

For  the  special  case  when  Bi  2  =  1 ,  o  =  a-\  =  o2  and  Equations 
(2-33a)  and  (2-33b)  can  be  simplified  to 

F(ara2)=2F(a)-2jf(u)F(u)du 

and  rearranging  finally  gives  the  joint  distribution  function  for  independent  and 
identical  case. (refer  to  Appendix  A-3): 

F(G1  ,o2)  =  2  F(o)  -  (F(o))2  (2-34) 

Equations  (2-33a)  and  (2-33b)  can  be  evaluated  for  specific 
models  such  as  the  Weibull  distribution  function. 

(a)  Weibull  Distribution  Model.   For  the  Weibull  model,  the  pdf  (f(x)) 
and  cumulative  density  function  (CDF:  F(x))  can  be  represented  by 


a-1 


f«HfXfJ  e*pHf 


a 


f| 


Substituting  Equation  (2-35)  and  (2-36)  into  Equation  (2-33a) 


(2-35) 
(2-36) 


F(ai'G2)  = 


o\ 


an      r 


1-exri-ly 


1-expUy 


a-\  -i 


^\m)M-m 


on  n 


1  -  expJ  - 


B12P 


du 


°2  a-1 


B 


a 


1  -  expi  - 


fB    v\ 

12 


a>  -i 


dv 
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By  integrating  and  rearranging  the  terms,  the  joint  distribution  function  of  the 
Weibull  model  under  bi-axial  stress  o^  and  o"2  can  De  simplified  to  (refer  to 
Appendix  A-4) 


a  ai 


^a1,a2)  =  1-exJ-[yj   -\-j 


(2-37) 


D.  MEAN  FOR  THE  JOINT  DISTRIBUTIONS 
1.   General  Case 

If  X-|  and  X2  are  two  continuous  random  variables  having  the  joint 
probability  density  function  f(o-|,02),  the  mean  or  expectation  of  X-|and  X2  can 
be  obtained  by: 

00         00 

^  =  E(°1)  =  1      J      01f(°V02)dCT1d02 


—   OO  —    OO 


(2-38a) 

^2  =  E(G2)  =  \      1      02f(01'°2)d01da2 


00  00 


■00     —  00 


(2-39a) 

Considering  the  stress,  the  bottom  limits  of  the  double  integration  are 
changed,  so  Equations  (2-38a)  and  (2-39a)  can  be  expanded  specifically  for 
the  stress  case  which  we  are  concerned  with 


00       00 


^1  =  E(0l)  =  J     Jaif(ai'°2)daidCT2 

v     y       0    0      v         *'  '  (2-38b) 


00       00 


^2  =  E(02)  =  J     /02f(<51'°2)d01d02  ,OOQ,v 

v      y       0    0       v  /  (2-39b) 

Similarly  if  three  random  variables  are  considered,  then  Equations 
(2-38b)  and  (2-39b)  can  be  expanded  by 

00       00      00 

h  =  J     J     J  <Jlf(<J1'a2'°6)dCT1dC2da6 

0    0    0      v  '  (2-40a) 

OO  OO  OO 

^2  =  J     \     J02f(ar02'°6)d01da2da6 

0    0    0    ^  v    '     ^     D/  *      D  (2-40b) 
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M-s  =  J     J     J  a6f(ai'a2'a6)daidCJ2da6 

o    o    o       v  '  (2-40c) 

2.   Specific  Case 

a.   Independent  Case 

If  the  random  variables  X-|  and  X2  are  independent  for  all  o-|  and  02 
f(a1>a2)  =  fi(a1)f2(a2) 
and  Equations  (2-38b)  and  (2-39b)  can  be  expanded  by 

00       00 

h-JB  Jn°1,l(Cl)f2(C2)d<rid02 

00       v    /     \     /  (2-41  a) 

OO  OO 

^2  =  1     l°2fl(ai)f2(°2)daidCJ2 

00         v     /     v      /  (2-42a) 

taking  the  integral  and  noting  that  F-j  (~)  =  F2(°°)  =  1,  then  Equations  (2-41  a) 
and  (2-42a)  can  be  simplified  to 

00 

^i  =  J°ifi(Gi)doi 

0        v     '  (2-41  b) 

00 

^2  =  J°2f2(a2)da2 

0         v     }  (2-42b) 

Similarly,  for  the  three  random  variable  case: 

00       00      00  00 

h"  J„   JQ  J0  01f,(CTl)f2(°2)f  6(°6)d01d02d(J6  =  J0O1fl(°l)dO1 

(2-43a) 

OO  OO  OO  OO 

^2  =  l0   J0  J0  a2fl(Gl)f2(a2)f6(G6)daida2d°6=l0a2f  2(a2)dCJ2 

(2-43b) 

00       00      00  00 

^6  =  J0   J0  J0  a6fl(ai)f2(a2)f6(°6)dO1da2dO6=l0°6f  6(a6)da6 

(2-43C) 
Observing  Equations  (2-41  b)  through  (2-43c),  it  can  be  concluded 
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that  the  mean  values  of  the  random  variables  X-j ,  X2,  and  Xq  are  dependent 
only  on  their  corresponding  probability  density  functions, 
b.  Independent  And  Identical  Case 

For  the  independent  and  identical  case,  the  joint  probability  function 
can  be  obtained  by 

f(a1)a2,a6)  =  f(ai)f(o2)f(a6) 
where  f(o)  =  i-\  (a)  =  f2(a)  =  1q(g) 

Applying  these  relations  to  Equations  (2-41  b)  through  (2-43c),  the 
mean  values  for  the  joint  distributions  can  be  simplified  for  the  independent  and 
identical  case. 


00 

=  Jof(ojd( 


U=U=U      =  JOT(CJ)ClG 

0  (2-44) 

(1 )  Weibull  Distribution  Model.    If   the  Weibull  probability  function  is 
considered,  Equation  (2-44)  can  be  expanded  as  follows 

oc-1 


-Uflf)  expKf)h 


letting    =  (o7p)a,  then 

[i  =  j    pyaexp(-y)dy 
0 

Using  the  definition  of  the  Gamma  function,  we  finally  get  the 
mean  function: 

*  =  K1  +  i)  (2-45) 

E.  JOINT  RELIABILITY 
1.   Concept  Of  Reliability 

The  common  notion  of  reliability  is  the  confidence  in  the  ability  of  a 
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device  to  perform  adequately  over  a  range  of  specified  conditions.  Such 
qualitative  measurement  of  confidence  and  adequacy  can  be  quantitatively 
calculated  by  probability.  In  the  formalism  of  probability,  the  conditions  over 
which  the  device  is  intended  to  perform  is  defined  as  the  random  variable  X: 
and  the  probability  over  the  range  of  conditions  is  defined  as  the  CDF  over  the 
range  that  the  random  variable  experiences  in  service.  That  is,  over  the  sample 
space  a  experienced  in  service,  the  probability  of  failure  is  defined  by  the  CDF 
over  the  sample  space: 
P(X  <  a)  =  F(g) 
Since  reliability  is  the  compliment  of  failure,  it  is  defined  by 

R(o)  =  1  -  F(o)  (2-46) 

2.   Joint  Reliability 

If  X-|  and  X2  are  two  random  variable,  then  the  joint  reliability  of  X-|  and 
X2  can  be  defined  by 

R(oi  ,02)  =  1  -  P(Xi  <oi  ,X2<a2)  (2-47) 

or  equivalents  for  the  three  random  variable  case 
R(o1,o2,o6)  =  1  -  P(X1<o1,X2<a2>X6<a6) 
For  the  discrete  case,  equation  (2-47)  can  be  expanded  by: 

usVa2  (2-48) 


F^ava2,a6)  =  1-  £    £    Xf(u,v,w) 


u<crv<a  „w<a  „ 

1  2  6 


If  the  continuous  case  is  considered,  the  joint  reliability  function  of  X-|  and 
X2  is  defined  by 


1   r      2 


Rfa1faJ)  =  1-f     J     f(u,v)dvdu 
v  J      J  (2-49) 


—  oo     — oo 


and  from  Equation  (2-49)  the  marginal  reliability  function  can  be  obtained: 
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1    °° 

RiCai)  =  1- J     ff  (u,v)dvdu 
A    V  J  (2-50a) 


—  oo     —  oo 


a 

2 


R?(°>?)  =  J[-!     Jf(u,v)dudv 
2^   2)  J     J  (2-50b) 


-  oo      —  oo 


If  the  three  random  variable  cases  are  considered,  Equation  (2-50a)  and 
(2-50b)  can  be  further  expanded: 

R/a  \  =  1  -  J  Jf  (u.v.w)dvdwdu 


—  oo     —  oo     — oo 


(2-51  a) 
R2(c2)  =  1-J  Jf  (u,v,w)dwdudv 

—  oo      —  oo     — oo 

o 

g       oo  oo 

R6(o6)  =  1-J     J     Jf  (u,v,w)dudvd 


oo  -oo  -oo  (2-51  b) 

a 


oo         oo 

'w 


R-l  (a-| ),  R2(o~2)>  ar|d  R6(°6)  are  a,so  equivalent  to  the  one  dimensional 
reliability  function  which  are  defined  by 

a 

R1(a.)  =  1-J   V(u)du 
A   V  _oo  1  (2-52a) 

R2(a2)  =  1-rf2(u)du  (252b) 

R6(a6)  =  1-J_\(u)du  (252c) 

a.   General  Case 

The  general  equations  for  reliability  can  now  be  applied  to  the 
calculation  of  the  mechanical  strength  reliability  of  composites  under  combined 
stresses  by  specifying  the  random  variables  as  Xj  and  the  realized  random 
variables  as  a\.  The  joint  reliability  under  combined  stress  g-\  and  02  can  be 
obtained  by 

F\ai'a2)  =  1  "  F(°i,a2)  =  1  "  P{X=X1  A  x20ver(ai'°2)} 


19 


where  X-|  and  X2  are  the  failure  strength  in  each  coordinate(see  Figure  2-2). 

Noting  that  the  joint  pdf  under  combined  stress  o-\  and  02  is 
represented  by  the  double  integration  of  f(a-| ,  02)  over  the  domains  S-|  and  S2 
and  that  the  double  integration  of  f(oi,  02)  over  the  domain  S3  (see  Figure  2-2) 
is  defined  as  the  joint  reliability  under  combined  stress  o^  and  02,  then  the 
following  relation  can  be  obtained: 


00       00 


FVc1fcr2)  =  1  -  F(ovg2)  =  J    J  f  (u.v)dudv 
v  7  v  J       °2  °i  (2-53) 

Comparing  Equation  (2-53)  to  Equation  (2-26),  Equation  (2-53)  can 

be  represented  in  another  form: 

f\ovo2)  =  Jl-j    J_iL_f(u,v)dvdu+  J       Jf(u,v)dudv 

0      b  °    R    v 

12  bi2V  (2-54) 

If  the  joint  reliability  under  three  combined  stress  is  considered, 

equation  (2-53)  can  be  expanded  by 


00       00       00 


F\aVa2»°6)  =  1  "  F(°1'a2'a6)  =  L    L   J0  f(U.V(W)dudvdw 

e     2     1  (2-55) 

b.   Specific  Case 

The  joint  reliability  of  bi-axial  stress  case  can  be  further  expanded  for 
the  specific  case  as  in  the  probability  distribution  function. 

(1)  Independent    Case.     If  the  random  variable  X-|   and  X2  are 
independent  for  all  o-j  and  02,  then  Equation  (2-53)  can  be  expanded: 


00       00 


f\°l.°2)=Uo,1(U),2(V)dUdV 

2     1  (2-56) 

After  integration,  Equation  (2-56)  can  be  simplified  to  (refer  to 
Appendix  A-5) 

R(a1,a2)  =  {1-F1(a1)}{1-F2(a2)}  (2-57) 

noting  that  R-|(o-|)  =  1-F-j  (a-| )  and  R2(o"2)  =  1"F2(°2) 
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R(o1lC2)  =  Ri(oi)R2(a2)  (2-58) 

For  the  three  random  variable  case,  Equations  (2-56),  (2-57),  and 
(2-58)  can  be  obtained  by 

oo        oo         oo 

F\ar°2'°6)  =  J0  J„  Jof,(u)yv)f6(w)dudVdw 

6        2        1 

=  (1-F1)(1-F2)(1-F6) 

=  R1(a1)R2(a2)R6(cT6)  (2-59) 

In  application  to  a  composite,  if  the  probability  of  failure  model  (Fj) 
is  known,  then  the  reliability  of  the  composite  can  be  calculated.  The  probability 
of  failure  model  can  be  inferred  from  the  physical  consideration  of  the  failure 
processes. 

When  a  filament  composite  is  loaded  along  the  fiber  direction, 
local  failure  begins  when  the  weakest  fiber  fails.  The  load  carried  by  the  broken 
fiber  is  transferred  to  the  neighboring  fibers.  Upon  additional  increase  in  load, 
additional  fibers  fail,  leading  to  the  increase  of  failure  sites  distributed  over  the 
composite.  The  higher  the  load,  the  higher  the  density  of  such  failure  sites  and 
the  higher  the  probability  of  clustering.  The  spatial  clustering  of  the  fiber  failure 
sites  leads  to  stress  concentration  and  ultimately  causes  the  catastrophic  failure 
of  the  composite.  Harlow  and  Phoenix  investigated  the  probabilistic  modeling 
of  the  above  sequential  failure  events  and  arrived  at  a  modified  weakest  link 
model  in  which  the  link  of  the  chain  is  a  bundle.  [Ref.  2]  For  a  limited  range  of 
the  random  variable,  it  can  be  approximated  by  the  two  parameter  Weibull 
model  [Ref.  3]: 


F^)  =  exp 


ot  1 


rUJ  I 


When  a  filament  composite  is  loaded  perpendicular  to  the  fiber 
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direction,  the  failure  mode  is  similar  to  that  described  by  classical  fracture 
mechanics.  The  exception  is  that  for  a  composite,  all  inherent  flaws  are  not 
randomly  oriented,  but  are  aligned  with  and  propagate  along  the  fiber.  In  such 
a  case,  the  largest  crack  dominates  and  it  is  effectively  a  weakest  link  process 
over  the  physical  volume.  On  these  physical  grounds  we  may  use  the  Weibull 
model  for  both  the  transverse  strength  and  shear  strength: 


F^a2)  =  exp< 


R(a6)  =  exp 


fo2\ 


2 


P 


K^2J 


^^ 


P 


V^6^ 


Equation  (2-57)  can  be  specified  for  these  Weibull  distribution 


functions: 


a  a 

1  2 


1 


Pi 


^ny 


Vh; 


F^ara2)  =  exp< 
equivalent^  for  the  three  random  variable  case: 
R(a1,a2,a6)=exp< 


(2-60) 


a  a  a 

1  2  6 

(cA\        fcn\         ^a„A 


k^j 


vP2; 


vhj 


(2-61) 


If  we  note    F(o-| ,  03)  =  1  -  R(a-) ,  02),  then  the  joint  distribution 
function  obtained  in  Equation  (2-32)  can  be  simplified  by 


F(oro2)  =  1-exp 


a  a  ^ 


1 


vPi; 


^2; 


(2-62) 


F(a1,a2,o6)  =  1  -exp 


"I* 


2  6 


Khj 


(2-63) 

(2)  Independent  And  Identical  Case.    If  the  random  variables  X-|  and 
X2  are  independent  for  all  c<\  and  02  and  the  probability  distribution  function  for 
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each  random  variable  is  the  same,  then  F(a)  =  F-((c-|)  =  F2(o"2).  Similarly  R(a)  = 
R-l  (a-|)  =  R2(<^2)  and  tne  Equation  (2-58)  can  be  simplified  to 

R(a1,o2)  =  R(a1)R(a2)  (2-64) 

As   a   specific   case,    if   the   Weibull   distribution    function    is 


considered,  Equation  (2-64)  can  be  further  expanded: 
^arc2)=exJ-(-£)    -(jf] 


(2-65) 

For  B-)  2  =  1.  then  a  =  c-\  =02  and  the  joint  reliability  function  for 
Weibull  model  is  obtained. 

R(a,o,  =  exp{-2(f)a}  ^ 

For  the  three  random  variable  case,  Equations  (2-64)  through 
(2-66)  can  be  expanded  by 

R(C1  ,a2,a6)  =  R(G1 )  R(a2)  R(a6)  (2-67) 


a  a  on 


R(a1,a2,a6)=exp|-^Tj    -  ^- j    -^ 


(2-68) 


on 
R(a,a,a)  =  exp<{-3[-£>) 

PV  J  (2-69) 
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III.  GEOMETRIC  REPRESENTATION  OF  JOINT  RELIABILITY  UNDER 

COMBINED  STRESS 

A.   BACKGROUND  FOR  GEOMETRIC  REPRESENTATION 

CDF  and  pdf  are  useful  functions  which  contain  all  the  relevant  information 
about  the  statistical  properties  of  a  random  variable,  in  our  case,  the  strength  of 
composite.  These  functions  are  required  in  reliability  analysis,  design, 
acceptance,  maintenance,  and  operational  logistics. 

In  applications,  in  order  to  map  out  the  statistical  failure  surface  even  for  the 
bivariate  case  (the  biaxial  combined  stress),  an  exceedingly  large  number  of 
experiments  are  required.  A  large  number  of  experiments  are  frequently 
impractical  due  to  economic  and  time  constrains  and  other  considerations. 
Therefore,  it  is  important  to  visualize  the  shape  of  the  failure  surface  in  order  to 
narrow  the  range  of  experiments  to  be  focused  on  the  critical  regions  in  the 
stress  domain. 

In  this  chapter  we  present  the  geometric  representations  for  the  reliability 
functions  for  the  independent  case  derived  in  Chapter  II.  The  graphical 
presentation  is  based  on  the  probability  plots  and  the  failure  surface 
representations,  which  were  also  investigated.  Examination  of  these  graphical 
representations  will  shed  light  on  the  appropriate  experiments  necessary  for  the 
identification  of  whether  the  failure  processes  are  independent  or  dependent. 
When  the  independence  is  established,  the  entire  failure  surface  (all 
permutations  and  combinations  of  combined  stress  case)  can  be  calculated 
from  the  uniaxial  strength  statistics. 
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B.   JOINT  RELIABILITY  FUNCTION 
1.   Two  Random  Variable  Case 

From  Equations  (2-56)  and  (2-66),  the  reliability  function  for  the  Weibull 
model  can  be  obtained  by 

VPJ  J  (3-1) 


R(a)  =  exp 


RELIABILITY 


E 

DC 


□     alpha=3 
♦      alpha=25 


sigma/beta 


Figure  3-1  :  Effects  Of  a  In  Reliability 
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If  we  consider  o7(3  as  a  normalized  variable,  then  the  reliability  function 
versus  o7{3  can  be  plotted  for  various  values  of  a.  Figure  3-1  shows  two  different 
graphes  for  large  and  small  values  of  a.  For  the  small  values  of  a,  the  graph 
shows  a  relatively  smooth  curve,  but  as  a  increases,  the  graph  approaches  to 
the  step  function  and  these  relations  can  be  further  expanded  for  the  joint 
reliability  problem.  If  the  independent  combined  stress  cases  are  considered, 
the  joint  reliability  of  the  composite  under  combined  stress  can  be  obtained  from 
Equation  (2-58). 

R(o1,02)  =  Ri(a1)R2(a2) 
Expanding  R(o-|,  o2) 


R(GrG2)  =  exp< 


a 


r1l 
s 

VH2  7 


(2-60) 


where 


R(°,)  =  exp'  - 


V0„ 


\R(<*2)  =  exp 


fa    \ 


8 

vK2y 


If  we  assume  that 

Pi  =Vi2P2 
and  noting  that 

o{  =  B12o-2 
then  the  joint  reliability  for  the  independent  case  can  be  represented  by 
R(a1,a2)  =  R1(a1)R2(a2) 


(3-2) 


=  exp 


a 


exp 


a 


B12^1 


(3-3) 


If  we  further  assume  Ci2=V-|2/Bi2.  then  Equation  (3-3)  can  be  simplified 


by 
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F^aro2)  =exp< 


a 


'exp 


-   C. 


2 


(3-4) 


From  Equation  (3-4),  the  reliability  functions  R-|(a-|)  and  R2(o"2)  can  De 
plotted  independently  in  two  dimensional  space  provided  that  a-|,  a2,  and  C-|2 
are  given.  Then  R-|(a-|)  and  R2(o"2)  are  combined  to  yield  the  joint  reliability 
function  of  R(a-|,a2)  as  shown  Figure  3-2.  Observing  the  Figure  3-2,  it  can  be 
noted  that  the  joint  reliability  is  affected  by  the  smaller  value  between  R-|(a-|) 
and  R2(^2)  f°r  a  Qiven  value  of  a-|/Pi .  that  is,  the  value  of  the  joint  reliability  is 
always  close  to  the  smaller  value  between  R-|  {a-\ )  and  R2(o~2)-  When  the  values 
of  both  R-|(a-|)  and  R2(o"2)  are  close  to  '1'  for  the  specific  range  of  g-|/Pi,  the 
value  of  the  joint  reliability  is  also  close  to  1,  but  as  either  R-|(a-|)  or  R2(o~2) 
decreases,  the  joint  reliability  also  decrease  depending  on  the  smaller  value 
between  R-|  and  R2-  We  can  also  observe  the  effect  of  C12  for  the  given  a-|  and 

cc2. 

As  shown  in  Figures  3-2a  and  3-2c,  the  R2(o"2)  shifts  to  the  right  from  the 
R-|(a-|)  as  C-|2  decreases  from  '1'  whereas  R2(o"2)  shifts  to  the  left  as  C-j  2 
increases  from  '1'  while  R-|(o"-|)  remains  constant.  So,  the  joint  reliability  follows 
the  smaller  function  between  R-|(ai)  and  R2(o"2)  and  appears  to  almost  overlap 
when  C-12  is  far  from  '1'.  As  the  value  of  C12  approaches  to  '1',  the  values  of 
R-|(g-|)  and  R2(o"2)  approach  each  other  and  when  C-|2=1-  Rl(°~l)  and  R2(o"2) 
cross  each  other  at  a-|/p-|=1.  The  joint  reliability  function  exists  to  the  left  of 
R-l  (a-| )  and  R2(<?2)  and  has  a  weaker  reliability  than  R-j  {g<\  )  or  R2(o"2)-  So  the 
joint  reliability  is  always  less  than  other  uni-axial  directional  reliabilities. 
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RELIABILITY  (A1=5,A2=3,C12=0.5) 
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Figure  3-2a  :  Reliability  Vs  oj/Pj  (C-j  2=0.5) 
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RELIABILITY  (A1=5,A2=3,C12=1) 
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Figure  3-2b  :  Reliability  Vs  aj/pj  (C-|2=1) 
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RELIABILITY  (A1=5,A2=3,C12=5) 
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Figure  3-2c  :  Reliability  Vs  g^  (C-|2=5) 
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Figure  3-3  shows  the  joint  reliability  in  three  dimensional  space  for 
different  values  of  C-|  2  and  for  fixed  a-j  and  0:2. 


A1=5,A2 


Figure  3-3  :  Reliability  In  3-D  Space  (aj=5,  a2=3) 
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Here  we  introduce  the  concept  of  linearized  failure  CDF  to  show  the 
effect  of  the  joint  failure  CDF  at  the  tail  area. 

To  get  the  linearized  failure  CDF,  denoted  by  F*(ct),  use  logarithmic 
algebra.  If  we  take  the  logarithm  of  R-|(o-|)  and  note  that  Ri(g-|)  =  1-F-|  (a-| ), 
then 


ln{Ri(°l)}  =  ln{1-Fi(°0}  = 


(3-5) 


Multiplying   both  side   by   (-1)  and  taking  the   logarithm  again,  the 
linearized  failure  CDF  can  then  be  defined  as 

Fi*(ci)-H-ln{Ri(ai)}J 


(3-6) 


Substituting  Equation  (3-5)  into  Equation  (3-6),  then 


F  *(<s  \  =  a.  In 
i  I    i)        i 


JiJ 


(3-7) 


By  the  same  way,  linearized  failure  CDF  corresponding  to  02  can  be 

obtained: 


F2*(a2)  =  a: 


KC12) 


+  ln 


K*1J 


(3-8) 


Comparing  Equation  (3-8)  to  (3-7),  it  can  be  noted  that  F2*(G2)  is 
vertically  shifted  and  rotated  when  compared  to  F-|*(o-|).  This  result  is 
represented  in  Figure  3-4. 

Similarly  F*(oi  ,02)  can  be  obtained  by 


F*(o1.a2)=lni 


( 


+ 


C   — 


(3-9) 


Comparing  Equation  (3-9)  to  (3-7)  and  (3-8),  it  can  be  noted  that 
F*(a-|  ,G2)  is  no  longer  a  linear  relation  in  terms  of  ln(oi  /pi ).   Figures  3-4a,  3-4b, 


32 


and  3-4c  show  the  linearized  failure  CDF  for  different  values  of  C-J2-  As  shown 
in  these  graphs,  F2*(a2)  is  shifted  upward  and  rotated  clockwise  as  C-\  2 
increases  and  this  causes  a  higher  probability  of  failure.  The  joint  linearized 
failure  CDF,  F*(ai,o2),  then  follows  the  higher  value  between  F-j  *(a-( )  and 
F2*(o~2)  and  does  not  show  the  linear  line,  especially  near  the  area  where 
F-|*(a1)  and  F2*(o"2)  intersect  each  other.  So  we  can  observe  that ,  as  the  value 
of  C12  increases  for  the  specified  cc-|  and  ot2,  F2*(a2)  shifts  upward  and  rotates 
clockwise  causing  the  intersection  between  F-j  *(a-| )  and  F2*(o"2)  to  go 
downward,  or  in  other  words  creating  a  higher  probability  of  failure. 

LINEARIZED  FAILURE  CDF 
(A1=5,A2=3,C12=0.5) 
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Figure  3-4a  :  Linearized  Failure  CDF  (C-|2=u-5) 
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LINEARIZED  FAILURE  CDF 
(A1=5,A2=3,C12=1) 
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Figure  3-4b  :  Linearized  Failure  CDF  (C-j  2=1 ) 
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LINEARIZED  FAILURE  CDF 
(A1=5,A2=3,C12=5) 


□  F1* 
♦  F2* 
+      F* 


In(sigma1/beta1) 


Figure  3-4c  :  Linearized  Failure  CDF  (C-|2=5) 
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As  a  special  case,  if  a=ai=a2,  then  the  joint  reliability  function  and 
Equations  (3-7)  through  (3-9)  can  be  simplified  by 


R(ai'G2)  =  eXP< 


(<5  ^ 


I    \\) 


0  +  Cn) 


F  *(G  \  =  aln 


(a^ 


F  *fc\  =  a< 


2*(*2) 


Jv 


a.  >\ 


F*(^2)=alnl^J  +  ln(1  +  Ci2) 


(3-10) 
(3-11) 
(3-12) 
(3-13) 


RELIABILITY  (A1=A2=5,C12=0.5) 
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Figure  3-5a  :  Reliability  Vs  ai/f3|(ai=a2=5) 
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RELIABILITY  (A1=A2=15,C12=0.5) 
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Figure  3-5b  :  Reliability  Vs  aj/pi  (a1=a2=15) 
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Comparing  Equations  (3-12)  and  (3-13)  to  Equation  (3-10),  it  can  be 
noted  that  F2*(o"2)  is  shifted  vertically  by  a  ln(v-|  2/B1 2)-   F*(a-|  ,02)  is  also  shifted 
vertically  by  ln{1  +(v-|  2/^1 2)}  regardless  of  values  of  C-|2-  Figures  3-5  though 
3-7  show  the  reliability  in  two  and  three  dimensional  space  and  linearized 
failure  CDF  in  two  dimensional  space  for  a  =  a-\  =  0:2- 


A1=A2=5 


Figure  3-6a  :  Joint  Reliability  In  3-D  (apao^) 
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Figure  3-6b  :  Joint  Reliability  In  3-D  (a1=a2=15) 
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As  shown  in  Figure  3-7,  each  failure  CDF  is  parallel  to  each  other  and 
the  joint  linearized  failure  CDF  almost  coincides  with  F-|*(o-|)  or  F2*(G2), 
depending  on  the  value  of  C12,  and  is  parallel  to  both  F-|*(o-|)  and  F2*(o2).  So 
in  the  physical  sense,  C.p  is  very  important  for  estimating  the  failure  and 
reliability.  The  joint  reliability  is  affected  by  the  weakest  reliability  function  and 
this  is  true  for  the  joint  linearized  failure  CDF. 
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(A1=A2=5,C12=0.5) 
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Figure  3-7a  :  Linearized  Failure  CDF  (aj=oc2=5) 
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LINEARIZED  FAILURE  CDF 
(A1=A2=15,C1 2=0.5) 
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Figure  3-7b:  Linearized  Failure  CDF  (a1=a2=15) 
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2.   Three  Random  Variable  Case 

If  the  three  random  variable  case  is  considered,  the  joint  reliability  under 
the  combined  stress  can  be  obtained  from  the  Equation  (2-79): 
R(a1 ,  a2,  a6)  =  R-|  (o^ )  R2(a2)  F^K) 
Expanding  this  equation  for  the  Weibull  model: 


f\Gvc2,c6)=  exp< 


a  a  a 

1  2  6 


vP2y 


vPey 


(2-61) 
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RELIABILITY  (A1=25,A2=5,A6=3) 
(C12=0.5,C16=1) 
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Figure  3-8a  :  Reliability  Vs  aj/Pj  (C-|  2=0-5,  C-|6=1) 
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RELIABILITY  (A1=25,A2=5,A6=3) 
(C12=1,C16=1) 
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Figure  3-8b  :  Reliability  Vs  ol/^l  (C-|  2=C-|  e=1) 
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If  we  assume  that 

P1=V12p2,(32  =  V16p6 
and  noting  that 

°1  =B12o2,  a-i  =B16o6 
where 

B     =^—  B     =  -J—  V     =-J—  V     =  — 

12       R      '16       R       '12      V       '16      V 


(3-14) 


(3-15) 


B     '12     V 

61  21 


'21  "61  '21  '61  (3-16) 

then  the  joint  reliability  for  the  independent  case  can  be  represented  by 


A1=25,A2  =  5,A6  =  3 


<j^ 


Figure  3-9  :  Joint  Reliability  In  3-D  (C-je=1) 
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R(ol,a2,a6)=R1  (a1)R2(a2)R6(CT6) 


=  exp 


=  exp 


1 


P. 


V.'  1/ 


exp 


exp 


v^2 ; 


C    -i 

S2p 


exp 


a 


exp 


a  >> 

C     -i- 


(3-17) 


where 


C     = 


V 


12 


12      B 


,C,   = 


V 


16 


12 


16       B 


16 


(3-18) 


Furthermore,  if  C12  and  ^16  are  given,  then  C26  can  be  calculated 

26        °6  02         °1CTeMl        B12V,6        C,6 


using  the  chain  rule: 
V 


C26"    | 


°2P, 


°2   ( 


1    PlP6"V12B16 


'26  ^  ^6  *       '    n    r6         "12  "16  12  (3-19) 

So,  if  a-j,oc2,oc6'  ^12«  ar|d  ^16  are  given,  then  the  joint  reliability 
function,  R(o-|  .o^.ag),  can  be  plotted  in  two  dimensional  space  with  respect  to 
o-|/(3-|.  Here  we  can  also  introduce  the  concept  of  linearized  failure  CDF  to  show 
the  effect  of  the  failure  CDF  at  the  tail  area.  To  get  the  linearized  failure  CDF, 
denoted  by  F*(o),  use  the  same  procedure  described  in  Equations  (3-5)  through 
(3-9),  then 


Pi 


vny 


F2>2)=aJln(C12)  +  ln 


r°,\ 


01 


\n  J 


V(G6)-a6{ln(C16)  +  ln(^} 


(3-7) 

(3-8) 

(3-20) 
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F*(o1Ia2,a6)  =  lrf 
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ra^\ 


p. 
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VC12R 
V         H  J 


<M  / 


V 


VC16B 


(3-21) 


So  it  can  be  noted  that  the  joint  linearized  failure  CDF,  F*(g-|  ,02>°6)  does 
not  show  a  linear  relation  in  terms  of  ln(a-|/pi)  except  when  a-|=a2=«6-  instead 
it  is  affected  by  F-j  *(a-| ),  F2*(o"2),  and  Fg*(og).  That  is  to  say,  F*(o-|,02,og)  is 
dominated  by  the  weakest  value  among  F-|*(a-|),  F2*(a2),  and  Fg*(ag)  so 
F*(o*-|  ,Q2,ag)  is  always  shown  on  the  top  of  each  of  graphs  as  shown  in  Figures 
3-1  Oa  and  3-1  Ob.  Figures  3-8a  through  3-8b  show  the  reliability  for  different 
values  of  C-|  2  when  a-| ,  0C2,  C-]  g  is  fixed.  These  figures  are  expanded  to  three 
dimensional  space  as  shown  in  Figure  3-9. 
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Figure  3-1 0a  :  Linearized  Failure  CDF  (C-|2=0-5.  C-j 6=1 ) 
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LINEARIZED  FAILURE  CDF  (C12=C16=1) 
(A1=25,A2=5,A6=3) 
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Figure  3-1  Ob  :  Linearized  Failure  CDF  (C-|2=Cl6=1) 
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As  a  special  case,  when  oc=a-|=a2=a5,  then  the  joint  reliability  and 
linearized  failure  CDF  functions  can  be  simplified  by 


a 

w  V  I     / 


V 


F6<oe>  =  0Yn(cie)  +  lltF 


\) 


(3-22) 
(3-23) 


F*ara2,a6)=alr(-^Uln(l-fC"2  +  C"6) 


VF1>  (3-24) 

Observing  Equations  (3-22)  through  (3-24),  it  can  be  noted  that  F-|*(a-|), 
F2*(a2)>  Fe*^)-  an<^  rT*(a-|  ,a2,cJ6)  are  parallel  each  other.  And  the  Joint 
linearized  failure  CDF,  F*(o-|  ,G2,c>q),  is  also  affected  by  the  largest  function 
among  F-|*,  F2*,  and  Fq*.  The  parameter  a  is  related  to  the  rotation  whereas 
C-|  2  and  C-|  6  are  related  to  the  vertical  shift  of  the  graph. 
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IV.  PROBABILISTIC  FAILURE  CONTOUR  UNDER  COMBINED  STRESS 

A.   PROBABILISTIC  FAILURE 

For  a  combined  stress  failure  mechanism,  when  one  stress  component  does 
not  affect  the  strength  of  another  component,  failure  is  mechanistically 
uncoupled,  as  shown  in  Figure  4-1.  When  one  stress  component  affects  the 
strength  of  another  component,  failure  is  then  mechanistically  coupled,  as 
shown  in  Figure  4-2. 


Figure  4-1  :  Mechanistically  Uncoupled  Mechanism 
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Figure  4-2  :  Mechanistically  Coupled  Mechanism 

If  the  probablistic  failure  is  considered,  the  mechanistically  uncoupled  case 
may  contribute  to  statistical  coupling  so  both  failure  mechanisms  are  not 
interacting  and  the  probablistic  failure  mechanism  is  distinguished  from  the 
deterministic  case.  In  order  words,  the  condition  for  the  failure  mechanism  is 
not  fixed,  it  depends  on  the  stochastic  combinations  of  the  intrinsic  strengths. 
So,  for  the  combined  stress  case,  uncoupled  independent  statistical  effects 
produce  phenomenologically  coupled  statistical  contours  and  the  mechanically 
coupled  dependent  statistical  effects  also  produce  phenomenologically  coupled 
statistical  contours.    In  the  following  subsection,  the  joint  probabilistic  failure 
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CDF  will  be  discussed  for  the  two  and  three  random  variable  cases  and  the 
effects  of  the  parameters  will  be  shown  graphically. 

B.  Two  Random  Variable  Case 

To  show  the  joint  failure  probability,  we  here  introduce  the  Weibull 
distribution  function.  From  Equation  (2-62),  the  joint  failure  distribution  function 
can  be  obtained  for  the  two  random  variable  case. 


F(a1,a2)  =  1-exp< 


'V 


a  a 


P, 


vny 


vP2y 


(2-62) 


From  the  above  equation,  (a^/p-j )  and  (02/P2)  can  De  computed  if  the  values 
of  F(o"i  ,02)  ancl  C12  are  9iven-  Figure  4-3  shows  the  joint  failure  CDF  in  three 
dimensional  space  and  Figure  4-4  shows  corresponding  joint  failure  contour  in 
two  dimensional  space  for  the  specific  value  of  F.  As  shown  in  Figures  4-3  and 
4-4,  the  shapes  of  the  failure  function  and  failure  contour  depend  on  the 
parameter  a,  which  is  material  dependent.  For  a  small  value  of  a,  the  failure 
contour  shows  a  smooth  curve  but  as  the  value  of  a  increases,  the  joint  failure 
contour  approaches  to  the  shape  of  a  rectangle.  This  phenomenon  is  the  same 
for  the  joint  reliability  contour  as  shown  in  Figure  4-5.  It  can  also  be  noted  that 
the  failure  contour  is  affected  by  a  large  value  of  a  between  a-j  and  0:2-  So  if 
anyone  of  the  values  of  a  is  larger  when  compared  to  another  value  of  a,  then 
the  failure  contour  approaches  to  the  shape  of  a  rectangle,  depending  on  the 
magnitude  of  a.  These  relations  are  shown  in  Figures  4-6  and  4-7.  Another 
result,  shown  in  these  failure  contour  graphs,  is  that  the  distance  between  the 
contour  line  is  small  when  a  is  large  compared  to  those  when  a  is  small.  So  the 
intersection  of  the  function  on  both  axis  depends  on  the  value  of  each  a. 
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A1=A2=5 


Figure  4-3a  :  Joint  Failure  CDF  In  3-D  (a1=a2=5) 
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A1=A2=25 


Figure  4-3b  :  Joint  Failure  CDF  In  3-D  (a1=a2=25) 


53 


1.50 


A1=A2=5 


1.00 


0.50 


0.1 


0.9 


0.00 


j i i i i i i i i i i_u i i_i i I  I  i     I     I     i     1   li     1     l     i     i     i     i     i     i     i 


0.00 


0.50 


1.00 


1.50 


Figure  4-4a  :  Joint  Failure  Contour  (a^=a2=5) 
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Figure  4-4b  :  Joint  Failure  Contour  (a1=a2=25) 
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Figure  4-5a  :  Joint  Reliability  Contour  (a1=oc2=5) 
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Figure  4-5b  :  Joint  Reliability  Contour  (a1=a2=25) 


57 


The  contour  line  of  F(o-|  ,02)  =0.5  exists  around  1.0  in  both  the  normalized  x 
and  y  axis  regardless  of  a,  but  as  a  increases,  the  F=0.5  contour  line 
approaches  to  1.0  in  both  normalized  axis. 
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Figure  4-6a  :  Joint  Failure  CDF  (aj=5,a2=3) 
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A1=25,A2  =  3 


Figure  4-6b  :  Joint  Failure  CDF  (a1=25,a2=3) 
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As  an  application,  if  we  compute  the  stress  of  a  specific  part  of  system  such 
as  a  rudder  or  ailerons  using  the  finite  element  method,  then  we  can  estimate 
the  probabilistic  failure  and  reliability  of  that  specific  part  of  the  system. 

The  parameters  a  and  p,  which  were  used  in  the  failure  CDF  and  reliability 
function,  are  material  dependent  so  these  values  should  be  determined  through 
experimentation  and  then  these  parameters  can  be  applied  for  that  specific 
material. 
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Figure  4-7a  :  Joint  Failure  Contour  (aj=5,a2=3) 
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Figure  4-7b  :  Joint  Failure  Contour  (a1=25,a2=3) 
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C.  THREE  RANDOM  VARIABLE  CASE 

If  the  three  random  variable  case  is  considered,  the  joint  failure  CDF  can  be 
developed  using  Equation  (2-63). 


F(Gva2,a6)  =  1  -exp 
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Figure  4-8a  :  Joint  Failure  CDF  (C-|  5=1 ) 
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A1=25,A2  =  5,A6  =  3 


Figure  4-8b  :  Joint  Failure  CDF  (C-|6=3) 
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From  the  above  equation,  normalized  axis  coordinates  can  be  computed  if 
the  values  of  F(o-j ,  o^?,  cig),  C-|  2,  and  C-|  6  are  given.  And  if  we  specify  the  value 
of  C-|6,  then  we  can  plot  the  joint  failure  CDF  in  three  dimensional  space  with 
respect  to  a-j/p-j  and  0-2/P2  as  shown  in  Figure  4-8.  As  the  value  of  C-|  6 
increases,  the  possibility  of  failure  increases,  that  is  to  say,  the  intersection  on 
the  normalized  x-axis  moves  to  the  left  because  oq/$q  term  is  embedded  in  the 
a-j/p-j  axis.  Figure  4-9  shows  the  corresponding  failure  contour  line  in  two 
dimensional  space  for  the  specific  value  of  F(o-| ,  o"2>  a6)- 
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Figure  4-9a  :  Joint  Failure  Contour  (C-j  6=1 ) 
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As  proven  graphically  in  the  three  random  variable  case,  the  probability  of 
failure  increases  as  the  number  of  random  variables  increase.  So  in  a  physical 
sense,  a  specific  composite  material  will  show  the  highest  reliability  when  the 
external  loads  are  applied  in  the  pricipal  axis  direction  because  there  are  only 
two  random  variable,  X-j  and  X2,  in  the  reliability  fuction.  But  when  the  external 
loads  are  applied  in  the  off-principal  axis  direction,  there  exist  some  shear  force 
random  variable,  Xq,  which  decreases  the  reliability  of  the  material. 
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Fig.  4-9b  :  Joint  Failure  Contour  (C-|2=3) 
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As  an  application,  if  the  stress  of  the  maneuvering  aircraft  is  computed  in  an 
arbitrary  direction,  then  the  stress  can  be  transformed  into  fiber,  shear,  and 
transverse  matrix  directions  which  can  then  be  used  to  predict  the  probability  of 
failure  at  that  specific  load.  Sample  calculations  were  made  and  the  results 
were  analyzed  with  respect  to  experimental  data  in  Appendix  B. 


66 


V.  CONCLUSIONS 

This  study  was  directed  towards  deriving  the  joint  reliability  and  joint  failure 
CDF  under  a  combined  stress  condition.  The  results  were  specialized  for  the 
Weibull  distribution  function  based  on  the  observations  that  the  composite 
material  failure  is  adequately  represented  by  the  weakest  link  model.  The  effect 
of  the  statistical  strength  parameters,  a  and  p,  which  are  material  dependent  are 
illustrated  using  two  and  three  dimensional  graphical  representations. 

To  analyze  the  reliability  of  the  composite,  the  inherent  statistical  strength 
parameters  in  the  composite's  principal  direction  need  to  be  experimentally 
measured.  Substitution  of  these  parameters  in  the  probabilistic  failure  criterion 
will  allow  for  the  estimation  of  the  reliability  of  the  composite  for  any  state  of 
stress. 

Comparison  of  joint  failure  distribution  under  the  restriction  of  independence 
to  experimental  data  suggests  that  mechanistic  coupling  of  the  failure 
mechanism  needs  to  be  included  in  future  extensions  of  the  formulation  and 
data  of  much  larger  number  of  samples  has  to  be  performed  at  critical  stress 
ratios  to  conclusively  examine  probabilistic  independence. 

Further  studies  may  be  extended  to  apply  these  reliability  and  failure  CDF  to 
the  specific  part  of  a  system  which  is  made  of  a  composite  material.  To  do  this, 
the  stress  of  the  specific  part  of  a  system  such  as  an  elevator  or  ailerons  should 
be  analyzed  using  such  as  finite  element  method.  These  equations  can  also  be 
used  in  a  repair  problem  which  requires  the  least  cost  and  most  effective 
method  to  analyze  the  combined  stress  in  a  specific  part  of  a  system. 
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So,  if  we  analyze  the  combined  stress,  the  appropriate  fiber  and  matrix  with 
proper  parameters  can  be  selected  to  make  the  composite  fit  to  that  specific 
part. 
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APPENDIX  A 
Chapter  II  Equations 

1 .   For  the  independent  case: 

P(X1+X2  :  X-i^a-i,  X2<o2)  =  P(X1<a1)  +  P{X2<02)  -  P(X-j^ai)  P(X2<cr2) 
=  Fi  (ai )  +  F2(o2)  -  F-|  (a-| )  F2(o2) 
If  the  Weibull  distribution  function  is  considered,  F-j  (a-j )  and  F2(o2)  can  be 
substituted  into  the  Weibull  distribution  function: 


P(X1  +  X2:X1<a1,X2<a2)  = 


1  -  exp* 


(g^ 


1 


1  -  exp* 


+ 


1  -  exp* 


fo,\ 


Pi 


1  -  exp 


V  H  J    )  ± 
By  expanding  and  simplifying 

P(X1  +  X2:X1  <  arX2  <  a2)  =  1  -  exp< 

which  is  the  same  result  as  the  Equation  (2-32) 
2.   For  the  independent  case: 

a 


U  JJ 


(o2\ 


\y2J 


a  a 


ro2\ 


a.   -\-\ 
2 


y^j 


vP2y 


F(ara2)  =  Jo  jj^j/u)  f2(v)dvldu  +  Jo  *lj~  ^(ujf  2(v)duldv 

I  12  J 

Taking  the  integrals  inside  the  parenthesis  first  and  substituting  limits  upon 
evaluation: 

F(<jv°2)=rfi(u>{F2(->-F2(^|du 
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+  Jo2f2(v){F1(-)-F1(B12v)}dv 
F1(a)  =  J°f  (u)du,F  (a)  =  Jf2(u)du 

o  o 

Noting  F-|  («>)  =  F2(°°)  =  1  from  the  definition  of  CDF  and  expanding 

o  a 

F(°r°2)  =  TVu>du- Jo  V^I^ 

o  a 

+  J    2f2(v)dv-J   'f    (v)F(B    v)dv 

o  o 

Taking  integrals  and  noting  that  F-|  (0)  =  F2(0)  =  0 

a 

F(°r°2)  =  F,(°0  +  F2l°2)  "  Jo  ,fi<u)F2(-§^)du 
"J   2,2(v>F,(B12v)dv 

o 

3.   From  Equation  (2-33a) 

Fcara2)  =  2F(a)-2jfa(u)F(u)du 

v  '  o 

By  integration  by  parts 


F(ara2)=2F(a)-2 


{F(a)}2-Jf(u)F(u)du 


noting 


J   f(u)F(u)du  =  0.5{F(o)}2 

o 

then  the  joint  distribution  function  for  independent  case  can  be  obtained  by 
F(d  ,o2)  =  2  F(a)  -  {F(a)}2 
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4.   By  expanding  and  integrating 


F(°i'a2)  = 


1  -  expi  - 
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on  -i 


1  -  expi  - 
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en  -i 
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_    2 


By  substituting  and  expanding 


F(ai'a2)  = 


1  -  expi  - 
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on  -i      r 
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on  n 


1  -  expi  - 
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+ 


expl- 


exp  - 
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f       Da         N 
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exp< 


a  on 

o\A       /   o\    \ 


pj    ^b12p 


-1 


exp' 
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12     2 
P       J 
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By  rearranging  terms  and  simplifying 


F(arG2)  =  1- 


Ba      ^ 

12 

X2+1v 


exp 


a  on 

vP~, 


lBi20y 


r  j ^ 

.B"  +1. 

V     12         J 


fon\ 


exp  - 


P 


vpy 


f      12    2 
P      J 


noting  a-|  =  B-|20~2,  tnen  tne  J°'nt  distribution  function  of  the  Weibull  model 
under  bi-axial  stress  a-j  and  02  can  be  simplified  by 
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a  on 


F(ara2)  =  1-exp|-^j 
5.   From  Equation  (2-50) 

f\0V°2)=CCf1<U),2<W»dudV 

2        1 
oo 

=L{Fi'"'-Ff-)}f2Mdv 

2 

noting  F-|(°°)  =  F2(°°)  =  1 

fVa1,a2)  =  {1-F1(a)}J<f2(v)dv 

^  ^  o 

integrating  then  we  finally  got  the  Equation  (2-57): 
R(a1,o2)  =  {1-F1(o)}{1-F2(a)} 
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APPENDIX  B 
Structural  Applications 

In  many  structural  applications  for  aircraft,  the  composite  laminae  are 
oriented  at  an  angle  to  the  tensile  load.  Take  for  example  the  skin  on  a  wing, 
where  the  filament  angle  is  oriented  primarily  along  the  direction  of  tension  due 
to  bending.  Transforming  the  stress  by  an  angle  9  to  an  off-axis  introduces  a 
shear  coupling  that  will  compensate  for  the  bending/rotational  coupling  of  a 
swept  wing.  Such  a  material  loading  relation  is  representative  of  that  indicated 
in  Figure  B-1.  We  will  explore  the  combined  stress  state  for  such  a 
configuration  and  the  structural  reliability  as  affected  by  the  off-axis  angle  0  . 
Since  we  are  addressing  a  spatially  homogeneous  state  of  stress,  both  stress 
analysis  and  the  substitution  of  the  combined  stress  in  the  probabilistic  failure 
criterion  can  be  carried  out  explicitly. 
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Figure  B-1  :  Transformation  Of  Coordinates 
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For  the  stress  analysis,  a  load  p  is  applied  to  the  composite  material  and  we 
denote  the  direction  of  p  to  be  x-|\  then  p  can  be  transformed  to  the  principal 
direction  of  the  composite;  xi  and  X2  where 

x-|  :  fiber  direction 

X2  :  transverse  matrix  direction 
and  the  angle  between  x-f  and  x-|-axis  can  be  defined  as  9  which  is  positive 
clockwise.  Using  the  tensor  relation,  the  stress  can  be  transformed  from  x-|'-X2' 
axis  to  x-|  -X2  axis  as  follows: 

°"jj  =  ami  anj  °mn'  (B-1) 

where 

i,  j,  m,  n  =  1 ,  2  for  2-dimension 
cn-p 

a22'  =  <J12'  =  <y2l'  =  ° 

ay  =  ^cos(x-|',  x-j)      cos(x-|',  X2)  \      =  /  cos(0)     -sin(0)  \ 

'  cos(x2',  x-j)      cos(x2',X2)  '         \  sin(6)      cos(0)  /  (B-2) 

Expanding  Equation  (B-1)  for  o11?  a^.  °"22' tnen 
°11  =am1  an1  °Ynn' 

=  311  an  an' +  an  a2i  ^12*  +  a21  ai2°"2l'  +  a21  a2l  <*22' 
=  a2i2  an' 
o11=pcos2(8)  (B-3) 

For  the  same  reason 

a12  =  -p  cos(0)  sin(0)  (B-4) 

o22  =  P  sin(0)  (B-5) 

To  calculate  the  reliability  associated  with  the  applied  load  p,  the  stress 

components  in  terms  of  p  (Equations  (B-3),  (B-4),  and  (B-5))  can  be  substituted 

into  Equation  (2-83)  to  find  the  joint  failure  CDF  for  the  specified  external  load  p. 
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We  note  that  Equation  (2-83)  is  specialized  for  the  case  that  the  effect  of  a^, 
aj2,  and  <5j2  on  tne  respective  failure  modes  are  independent.  Substitution 
yields: 


F(p,0)  =  1-exp 


pcos^(O) 
P, 


a  a 


a 


psirr(9) 

To 


fpcos(0)sin(0n 


P, 


(B-6) 


where 

al=all'02  =  a22'a6  =  a12 
Comparing  Equation  (B-6)  to  (2-83)  and  (3-18),  it  can  be  noted  that 
B12  =  cot2(0),  B16  =  -cot(0) 
Ct  2  =  V-|  2  tan2(8),  d 6  =  -Vi  6  tan(0) 
This  joint  failure  CDF  in  terms  of  p  and  0  can  be  plotted  in  terms  of  p  for  an 
specified  0.  Similarly,  from  the  Equation  (B-6),  the  linearized  failure  CDF  can  be 
obtained  by 


F"  (p,  6)  =  I  n 


a 


a 


pcos  (0) 


psin2(0) 

Po 


pcos(0)sin(0) 


(B-7) 


In   order  to    investigate   the   dependency   of  the   failure    mechanism, 
experimental  data  from  Sun  and  Yamada  [Ref.  4]  is  examined.     From  the 
experiment,  the  following  parameters  were  obtained  for  fiber,  transverse,  and 
shear  force  using  uni-axial  tests. 
OCJ  =20.5,(3!  =  127000 
oc2  =  5.66,  p2  =  1070 
a6  =  8.96,  P6  =  3078 
Substituting  these  strength  parameters  into  Equations  (B-6)  and  (B-7),  the 
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probability  of  failure  for  load  p  applied  to  different  angles  can  be  examined. 
Figures  (B-2a)  through  (B-6b)  compare  the  joint  failure  CDF  to  the  experimental 
data  for  different  angles.  These  data  will  be  examined  individually  for  each  of 
the  five  angles  where  experimental  data  is  available. 
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Figure  B-2a  :  Comparison  Of  Experimental  Data  with 
Joint  Failure  CDF  For  Independent  Case  (15  Deg.) 
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At  8  =  15°,  even  when  the  load  is  closely  aligned  with  the  fiber  angle,  Figure 
B-2a  shows  that  the  fiber  strength  statistics  has  no  effect  on  the  combined 
failure  probability  by  noting  that  the  fiber  failure  curve  is  off-scale  for  the  high  p 
region.  That  is,  the  failure  is  always  dominated  by  a  combination  of  matrix 
transfer  strength  and  the  matrix  shear  strength.  The  interaction  of  transverse 
and  shear  can  be  better  observed  in  the  linearized  Weibull  plot,  Figure  B-2b, 
from  which  it  can  be  seen  that  in  the  lower  tail  (low  p  region)  the  transverse 
strength  dominates  and  in  the  upper  tail  (high  p  region)  the  shear  strength 
dominates. 
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Figure  B-2b  :  Comparison  Of  Experimental  Data  With  Joint 
Failure  CDF  For  Independent  Case  (Linearized,  15  Deg.) 
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While  the  experimental  data  are  not  close  to  the  predicted  joint  failure  CDF 
under  the  simplification  of  independence,  its  significance  can  not  be  evaluated 
because  of  the  small  sample  size.  With  the  eight  samples  in  the  current  case, 
neither  the  upper  tail  nor  the  lower  tail  can  be  observed.thus  the  failure 
probability  coupling  remains  unresolved.  The  location  discrepancy  perhaps 
has  more  significance.  It  suggests  that  combined  transverse  and  shear 
stresses  have  weakened  the  composites  and  that  the  failure  mechanistic 
coupling  needs  to  be  treated. 
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Figure  B-3a  :  Comparison  Of  Experimental  Data  With 
Joint  Failure  CDF  For  Independent  Case  (20  Deg.) 
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At  0  =  20°,  Figures  B-3a  and  3b  show  that  even  with  a  slight  angle  change  of 
5°,  the  effect  of  the  shear  component  diminishes  rapidly.  The  transition  region 
has  shifted  higher  on  the  upper  tail.  Physically,  it  means  that  unless  a  large 
number  of  samples  are  tested,  only  the  strongest  of  the  samples  will  fail  in 
shear,  and  the  remaining  will  fail  by  transverse  stress.  For  this  reason,  and 
again  because  of  the  small  number  of  samples,  the  shape  indicated  by  the 
experimental  data  has  no  significance.  The  location  difference  between  the 
data  and  prediction  suggests  that  mechanistic  coupling  needs  to  be  included. 
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Figure  3-3b  :  Comparison  Of  Experimental  Data  WithJoint 
Failure  CDF  For  Independent  Case  (Linearized,  20  Deg) 
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At  0  =  30°,  Figures  B-4a  and  4b  show  that  the  effect  of  shear  is  shifted  to 
even  a  higher  upper  tail.  All  comments  on  0  =  20°  data  apply  to  0  =  30°  as  well. 
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Figure  B-4a  :  Comparison  Of  Experimental  Data  With 
Joint  Failure  CDF  For  Independent  Case  (30  Deg.) 
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Figure  B-4b  :  Comparison  Of  Experimental  Data  With  Joint 
Failure  CDF  For  Independent  Case  (Linearized,  30  Deg.) 
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At  0  =  45°,Figures  B-5a  and  5b  show  that  the  effect  of  shear  is  practically 
nonexistent;  in  fact,  it  can  not  be  observed  from  the  range  of  the  scale 
presented.  What  is  significant  is  the  location  of  data  is  below  the  location  of  the 
joint  CDF  suggesting  mechanistic  coupling. 
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Figure  B-5a  :  Comparison  Of  Experimental  Data  With  Joint 
Failure  CDF  For  Independent  Case  (Linearized,  45  Deg.) 
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Figure  B-5b  :  Comparison  Of  Experimental  Data  With  Joint 
Failure  CDF  For  Independent  Case  (Linearized,  45  Deg.) 
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At  0  =  60°, Figures  B-5a  and  5b  show  that  the  effect  of  shear  is  no  longer 
present  and  that  the  fit  of  the  data  is  much  improved.  The  latter  observation 
further  substantiates  the  existence  of  mechanistic  coupling. 
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Figure  B-6a  :  Comparison  Of  Experimental  Data  With 
Joint  Failure  CDF  For  Independent  Case  (60  Deg.) 
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Figure  B-6b  :  Comparison  Of  Experimental  Data  With  Joint 
Failure  CDF  For  Independent  Case  (Linearized,  60  Deg.) 
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From  the  experimental  data  available  in  the  open  literature,  it  was  observed 
that  a  much  larger  number  of  samples  are  required  to  identify  the  probabilistic 
coupling  and  dependency  of  strength.  What  is  evident  is  that  mechanistic 
coupling  needs  to  be  included  as  an  extension  to  the  investigation  herein. 
Furthermore,  experimental  design  using  the  probability  failure  criterion  is 
mandatory  to  optimize  experimental  testing. 
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